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Abstract: We present the results of a long M87 monitoring campaign in very high energy 7-rays with the MAGIC-I 
Cherenkov telescope. A total of 150 hours of data was gathered between 2005 and 2007. No flaring activity was found 
during that time. Never the less, we have found an apparently steady and weak signal at the level of 7a. We present the 
spectrum between 100 GeV and 2 TeV, which is consistent with a simple power law with a spectral index —2.21 ±0.21 
and a flux normalization (at 300 GeV) of 5.4 ± 1.1 x 10 ~^ xeVsm^ • complements well with the previously published 
Fermi spectrum, covering an energy range of four orders of magnitude without apparent change in the spectral index. 
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1 Introduction 

M87 is a giant elliptical radio galaxy situated in the Virgo 
cluster at a distance of 16.7 Mpc |T|. The centre of the 
galaxy is a SMBH with a mass of (6.4 ± 0.5) x lO^M© 
El . The prominent jet of M87 (first discovered in the op- 
tical band O) is inclined by 10° — 45° from our line of 
sight The jet is not homogeneous but shows several 

distinctive regions. The compact nucleus has a dimension 
of several hundreds of Schwarzschild radii. Throughout 
the inner jet several brighter spots, so called "knots", are 
visible, often characterized by superluminal motion with 
apparent speeds reaching ~ 6c |4|. The jet continues to 
expand into the inter-cluster medium, forming giant lobes 
with a length of 80 kpc. A closer inspection also indicates 
the existence of a counter-jet. The large scale jet primarily 
emits from radio to X-rays through the synchrotron mech- 
anism. The two most promising emission regions at all 
wavelengths (including very high energy 7-rays) are the 
nucleus and the knot closest to the nucleus (60 pc), known 
as HST-1. 

The first indication that M87 is a source of VHE 7-radiation 
was reported by the HEGRA collaboration |6|, which was 
confirmed by the H.E.S.S collaboration with data collected 
between 2003 and 2006 | 7 1. This contribution will focus on 
long-term monitoring observations of M87 with MAGIC 
during a low emission state in 2005-2007. 



2 Observations and Data Analysis 

MAGIC started to observe M87 regularly in 2005 and con- 
tinued to monitor it during several observational cycles un- 
til today (see e.g. the recent announcement of an enhanced 
emission state in February 2010 |8|). This paper focuses 
on data taken between March 2005 and June 2007 with the 
stand-alone MAGIC-I telescope. The MAGIC-I telescope 
underwent several significant upgrades over this period. 
Specifically these upgrades include: The optical PSF of the 
17m mirror dish improvement. The readout chain has been 
significantly enhanced during the observation time. The 
data discussed in this paper include three different readout 
stages: 1) 300MS/s readout without splitters, 2) 300MS/s 
readout including splitters and 3)2 GS/s readout with split- 
ters. After the installation of the 2 GS/s Mux-FADC system 
a detailed study of the timing properties of the showers lead 
to a significant improvement of the sensitivity of the tele- 
scope |9|. 

The analysis technique used throughout this paper is the 
same one performed in ifTOl , where an excellent agree- 
ment between the published Crab Nebula spectrum and the 
stacked analysis was found, confirming the stability of the 
analysis over large time scales. 

MAGIC observations of M87 were performed in two dif- 
ferent modes: ON-OFF and wobble. In the ON-OFF mode 
the telescope is pointed directly at the source during the 
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Table 1 : Observation periods 



ON runs. The background is estimated with dedicated OFF 
data. In the so-called wobble mode |11| the telescope is 
not pointed directly at the source position but 0.4° away. 
Accordingly wobble data does not require dedicated OFF 
runs. Only the 2005 observations of M87 were recorded 
in ON-OFF mode while the majority of the data taken in 
the later years was taken in wobble mode. Throughout this 
paper we have applied the same timing based analysis to 
the entire dataset, which is described in detail in |[T2ll . Each 
period with different telescope conditions was first anal- 
ysed separately with matching OFF-data (when needed) 
and Monte Carlo simulations from the corresponding pe- 
riod and finally combined after the 7-hadron separation. 

The entire observation campaign consists of 154.1 h out of 
which 128.6 h survive the quality selection cuts. This se- 
lection ensures that only data with good weather conditions 
and without any technical problems is used in the analysis. 
Table [T] gives an overview how this data is distributed over 
the entire campaign. 

3 MAGIC results 

The 0'^ distribution of the entire observation campaign is 
shown in Fig. [T] An excess of 862 events over 12454 back- 
ground events is apparent, corresponding to a significance 
of 7.0 cr (calculated with formula 17 of lfT3l ). 

As the light curve in Fig. [3]shows, the excess over the entire 
observation campaign is compatible with a fit to a constant 
flux with a reduced of 0.51. Additionally we analysed 
each observation night independently, however none of the 
nights yielded an excess >3cr, which is consistent with ran- 
dom fluctuations. Fig. [2] shows the cumulative excess dis- 
tribution over time. The excess grows linearly, which is 
consistent with a constant emission. We thus conclude that 
no significant flare occurred during our observation cam- 
paign. Note that the data in 2005 are compatible within er- 
rors with the high flux state observed by H.E.S.S. |7|. The 
reason is mainly a too short MAGIC observation time (and 
accordingly large error on the source flux) as well as the 
short time variability of the source (flares on a time scale 
of two days have been reported by the H.E.S.S. Collabora- 
tion). 
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Figure 1: distribution of the combined 128.6 h MAGIC 
observations between 2005 and 2007. The grey shaded 
area below the red crosses is the OFF data-sample, while 
the black crosses corresponds to the ON data. The excess 
is point-like compared to the PSF with a significance of 7 
standard deviations (calculated with formula 17 of 1 131 ). 
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Figure 2: Evolution of excess events from M87 over time. 
Only wobble data are shown in this figure in order to reduce 
systematic uncertainties in the ON-OFF subtraction. 



The MAGIC spectrum is consistent with a simple power 
law with a spectral index -2.21 ±0.21 and a flux normaliza- 
tion (at 300 GeV) of 5.4 ± 1.1 x 10" %^^^^^, . Both the 
Fermi spectrum (taken from |14|) and the MAGIC spec- 
trum are shown in Fig. (U covering an energy range of four 
orders of magnitude without apparent change in the spec- 
tral index. 

A detailed discussion of the potential systematic errors can 
be found in lfT5l . As a summary the systematic energy 
scale error is about 16%, the error of the flux normaliza- 
tion is estimated to be 1 1 % and the systematic slope error 
is ±0.2. Using Monte Carlo simulations for each observa- 
tional epoch we find that the error of the flux normalization 
should be < 10% and the error of the slope should be <0.03 
(note that this applies only to hard spectra with a spectral 
index < 2.5 as is the case for the data discussed here). 



4 Modelling the SED 

We model the SED using the structured jet model of (161, 
already applied to M87 in y/7J. The model assumes that the 
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Figure 3: Light curve of the integral 7-ray flux above 
100 GeV during the 2005-2007 MAGIC observation pe- 
riod. The dashed blue line corresponds to the fit re- 
sult of a linear function to the data points, with a re- 
duced x^/d.o.f. = 5.66/11 = 0.51 and a mean flux of 
5-06x10-12^^. 
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Figure 4: Shown above is the combined Fermi-MAGIC 
differential energy spectrum over four orders of magni- 
tude in energy starting from lOOMeV up until 1.5 TeV. The 
dashed line corresponds to the extrapolation of the fit to the 
MAGIC data points into the Fermi energy range. Note that 
the two observations are not contemporaneous. 

jet has an inner fast core (spine), with bulk Lorentz factor 
Fs, surrounded by a slower layer, with bulk Lorentz fac- 
tor Fl. In both regions, relativistic electrons emit through 
synchrotron and inverse Compton mechanisms. 

The model is fully specified by the following parameters: 
i) the spine is assumed to be a cylinder of radius R, height 
Hs (as measured in the spine frame) and in motion with 
bulk Lorentz factor Fs; ii) the layer is modelled as a hol- 
low cylinder with internal radius R, external radius R2, 
height i^L (as measured in the frame of the layer) and bulk 
Lorentz factor Fl. Each region contains tangled magnetic 
field with intensity Bs, Bi, and it is filled by relativistic 
electrons assumed to follow a (purely phenomenological) 
smoothed broken power-law distribution extending from 
7min to 7max and with indices ni, n2 below and above 
the break at 7b. The normalisation of this distribution is 
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Figure 5: SED of the core of M87 (green open squares) 
together with the MAGIC data points (black filled circles). 
Refer to the text for details. 

calculated assuming that the system produces an assumed 
(bolometric) synchrotron luminosity Lgyn (as measured in 
the local frame), which is an input parameter of the model. 
We assume that Hi, > Hs. The seed photons for the IC 
scattering are not only those produced locally in the spine 
(layer), but we also consider the photons produced in the 
layer (spine). 

The result of the modelling is reported in Fig. [51 where 
we show the SED of the emission produced by the spine 
(red) and the layer (blue) and their sum (black). Green 
open symbols report historical data, black filled circles 
show the MAGIC spectrum. For comparison we also re- 
port the H.E.S.S. spectra taken in 2004 and 2005 (blue and 
red circles, respectively, from JTl) and the 2007 VERITAS 
spectrum (cyan, from |18|). Clearly the 2005-2007 spec- 
trum was very similar to that measured by H.E.S.S. in the 
"low" 2004 state and by VERITAS. The green "bow-tie" 
reports the nuclear X-ray spectrum measured by Chandra 
in 2000 [19]. The black "bow-tie" show the average X-ray 
spectrum during the period covered by the MAGIC obser- 
vations, obtained by scaling the 2000 spectrum by a factor 
of 4, as inferred in from the lightcurve reported by [20|. 
During the MAGIC observing periods there are no mea- 
surements in the GeV energy band. However, for compar- 
ison, we report the Fermi/LAI spectrum obtained integrat- 
ing over the first 10 months of all-sky survey data |[T4ll with 
a black bow-tie. The green open square in the high energy 
7-ray range corresponds to the EGRET upper limit. The 
SED shows two pronounced bumps, one peaking in the 
IR band, one extending from MeV to TeV energies. The 
adopted parameters are reported in Table [21 

A feature of the realization of the "spine-layer" model for 
M87 is that the IC emission from the layer, accounting 
for the observed VHE emission, is partly absorbed through 
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Table 2: Input parameters of the models for the layer and the spine shown in Fig. [51 All quantities (except the bulk Lorentz 
factors r and the viewing angle 0) are measured in the rest frame of the emitting plasma. The external radius of the layer 
is fixed to the value R2 = 1.2 x R. 



the interaction with the optical-IR photons produced in the 
spine. This leads to a relatively soft spectrum, suitable to 
reproduce the low level spectrum measured by MAGIC, 
but difficult to reconcile with the hard spectrum recorded 
during high state by H.E.S.S. A possible way out of this 
problem is to enlarge the emission regions. However, the 
increase of the source size is limited by the observed short 
variability time- scales. 

It appears that the power associated to the magnetic field is 
largely dominant over that related to relativistic electrons 
and protons (assuming a composition of 1 cold proton per 
emitting electron, see e.g. |21 1) both in the jet and in the 
layer. In any case the jet power required to reproduce the 
SED, strongly dominated by the spine, appears relatively 
modest, Pjet = 1.5 x 10^^ erg/s. 

5 Summary 

MAGIC has detected weak and steady VHE 7-ray emission 
from M87 between 2005 and 2007. Our measurements are 
compatible with previously reported low states by VERI- 
TAS and H.E.S.S., which suggests that the observed emis- 
sion level and spectral characteristics are extremely stable 
and may refer to the "ground- state" emission of M87. This 
is also supported by the extrapolation of the MAGIC spec- 
trum into the Fermi energy range, which is in perfect agree- 
ment despite the fact that both observations are not contem- 
poraneous. 

We were able to describe this emission with a structured 
jet model, that separates the jet into a spine and an outer 
layer. An interesting feature of the model is the existence 
of a transition region between the emission from the spine 
and that from the layer around 40 GeV. Currently no mea- 
surements of M87 are available in this energy range. This 
is due to the small collection area of the Fermi satellite and 
the insufficient sensitivity of Cherenkov telescopes in the 
low energy domain. Since fall 2009 MAGIC consists of 
two 17 m diameter telescopes, which observe in a stereo- 
scopic mode. This upgrade has especially improved the in- 
strument's capabilities at energies below 100 GeV |[22l. It 
is thus possible that a future, deep observation of M87 with 
the stereoscopic MAGIC system will reveal a feature in the 
otherwise smooth power law spectrum of M87, confirming 
or denying the validity of the structured jet model. 
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